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ABSTRACT 



Aims. We revisit the scaling relations and star-forming histories of local elliptical galaxies using a novel selection method applied to 
the Sloan Digital Sky Survey DR7. 

Methods. We combine two probability-based automated spectroscopic and morphological classifications of ~ 600000 galaxies with 
z < 0.25 to isolate true elliptical galaxies. Our sample selection method does not introduce artificial cuts in the parameters describing 
the galaxy but instead it associates to every object a weight measuring the probability of being in a given spectro-morphological class. 
Thus the sample minimizes the selection biases. 

Results. We show that morphologically defined ellipticals are basically distributed in 3 spectral classes, which dominate at different 
stellar masses. The bulk of the population (~ 50%) is formed by a well defined class of galaxies with old stellar populations that 
formed their stars at very early epochs in a short episode of star formation. They dominate the scaling relations of elliptical galaxies 
known from previous works and represent the canonical elliptical class. At the low mass end, we find a population of slightly larger 
ellipticals, with smaller velocity dispersions at fixed stellar mass, which seem to have experienced a more recent episode of star 
formation probably triggered by gas-rich minor mergers. The high mass end tends to be dominated by a third spectral class, slightly 
more metal rich and with more efficient stellar formation than the reference class. This third class contributes to the curvature of 
the mass-size relation at high masses reported in previous works. Our method is therefore able to isolate typical spectra of elliptical 
galaxies following different evolutive pathways. 

Key words. Galaxies: fundamental parameters, Galaxies: evolution, Galaxies formation, Galaxies: elliptical and lenticular, cD 



1. Introduction 

More than 50% of the stellar mass in today's universe is 
in early-type galaxies (see Fukugita et al. 1998| ). Understanding 
their formation and evolution is therefore one of the key ques- 
tions in modern structure formation theories. 

Basically, there is an open debate on whether actual elliptical 
galaxies are the end products of a long hi story of major and/or 
minor dry/wet galaxy mergers (see e.g. IDe Lucia et al. 20061 
hierarchical scenario) or if they are passively evolving galax- 
ies formed at high redshift in a short time- scale with al- 
most no accretion of external material since their early 
formation (monolithic scenario; see e.g. [Kodama et al~1 998, 
Ivan Dokkum & Ellis 20031 ICarretero et al. 20071) . From the ob- 
servational point of view, the scaling relations are a powerful 
tool to disentangle between these two evolutionary paths. 

The passively evolving scenario or monolithic collapse 
model ( |Partridge & Peebles 1967| ILarson 19751 ) is supported 
by the fact that early-type galaxies form, at first sight, an 
homogeneous family of objects with very regular morphol- 
ogy (they do not show strong spiral arms or asymmetries 
like other types of galaxies; e.g. Abraham et al. 1996). In 
addition, their surface brightness distribution is concentrated 
towards the center, and can be well fitted by a Sersic surface 



brightness profile (ICaon et al. 19931 jTrujillo et al. 2001 
Gra ham& Guzman 20031 |Aguerri et al. 2004 

Gutierrez et al. 20QH 

|Kormendy et al. 2009) . This regular morphology is also 
reflected in the color (i.e., stellar populations). Early-type 
galaxies form indeed a well defined family of galaxies in 
the color-magnitude diagram on the so-called red sequence 
(IKodama & Arimoto 19971 IBell et al. 20041 |Balogh et al. 2004| ) 
which suggests that they are formed by old stellar popula- 
tions of similar age with almost no star for mation activity at 
present. This holds at le ast up to redshift 1 (IMei et al. 20091) 
and even up to redshift 2 ( Andreon & Huertas-Company 201 1| ). 
From the dynamical point of view, they are dominated 
by random motions with a high degree of virializa- 
tion (see e.g. IKronawitter et al. 20001 IGerhard et al. 20011 
IThomas et al. 2007T) . They consequently follow tight 
scaling relations defined by their size, velocity dis- 
persion, surface brightness, and luminosity. The most 
comm only studied r el ations are: the lu minosity/mass- 
size (IShen et al. 2"003l IMcIntosh et al. 2005)) . the Faber- 
Jackson relation (IFaber & Jackson 19761 ), the Kormendy 
relation (|Kormendy 1977 ), and the fund amental plane 



(Djorgovski & Davis 1987 



Dressier et al. 19871) . 
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Detailed studies on scaling relations of early-type galaxies 
have shown that despite of the apparent homogeneity, there are 
important differences among them. Some works have in fact 
found that the scaling relations and/or structural parameters 
of early-type galaxies depend on different galaxy properties, 
such as: luminosity dBender et al. 19921 [C aon et a l. 19931 
Desroches et al. 20071 |Nigoche-Netro et al. 2008 



Aguerri & Gonzalez- Garcia 2009|Nigoche-Netro et al. 2009 



Nig oche-Netro et al. 2010|Bernardi et al. 2010b" 

stellar populations dForbes et al. 1 998, 

ITerlevich & Forbes 2002IGraves et al. 20091 

Graves & Faber 2010), or environme nt ( [Aguerri et al. 2 004, 
IGutierrez et al. 200"4l IFritz et al. 2005b . Moreover, a strong 
size evolution is observed in early-type galaxies from high 
to low redshift (IDaddi et al. 20051 jTrujillo et al. 2006[ ). These 
differences can be only understood if galaxies with different 
luminosities, stellar populations or located in different environ- 
ments have experienced different formation histories. According 
to these properties, some authors have seen a dichotomy in 
the early-type population, and divided them in two different 
families. One of the families would be formed by the low lumi- 
nous early-type galaxies. They show fast rotation, are isotropic 
oblate spheroids with discy isophotes, and show flat surface 
brightness profiles in their central regions (they are core-less 
objects). These galaxies are similar to bulges of spiral galaxies. 
In contrast, the more luminous or giant early-type galaxies 
are non-rotating, anisotropic and triaxial systems. In addition, 
they show cuspy cores and boxy isophotes. These galaxies are 
also older and more a-enhanced than the less luminous coun- 
terparts ([Kormendy & Illingworth 1982[ INieto et al. 1991 



IBenderetal. 19921 
IFaber et al. 19971 
lEmsellem et al. 20071 



Kormendy & Bender 1996 
ILauer et al. 20051 ILauer et al.2007 



|Kormendy et al. 2009[ ). 
lenged by other 



|Cappellari et al. 2007 



This dichotomy has been chal- 
authors ([Jerjen & Binggeli 1997 



IGraham& Guzman 20031 IGavazzi et al. 2005b . 

The differences reported above in early-type galaxies 
could be explained if their formation history is some- 
how more complex than what is suggested by the mono- 
lithic scenario. Since the pioneering work by Toomre & 
Toomre (1972), it is thought that early-type galaxies can 
be formed by mergers of disk galaxies. Detailed simula- 
tions report that two equal disks do no t survive after a 
merger process (|Barnes& Hernquist 1991 
Barnes & Hernquist 1996| 



Hernquist 1993 
|Mihos & Hernquist 1996 



Pi Matteo et al. 2005, Naab et al. 2006b, Robe rtson et al. 2006 



Hopkins et al. 2009). The end product of such mergers is 
a galaxy with early-type morphology. These numerical 
simulations have also shown that the differences observed 
in the properties of early-type galaxies can be explained 
through the different physical processes taking place in such 
kind of mergers. This involves gas dissipation, star forma- 
tion, and super-massive black hole feedback dBarnes 19921 
Barnes & Hernquist 1996[ |Mihos & Hernquist 1996 1 



Di Matteo et al. 2005, Robertson et al. 2006, Naab et al. 2006a 
Naa b et al. 2006bl INaab et al. 20071 |Hopkins et al. 2008| ). In 
particular, the gas fraction in the galaxy progenitors could 
play a crucial role in shaping the final merger remnants and 
fix their scaling relations (see, e.g., IRobertson et al. 2 006, 
|Hopkins et al. 2008 ). 

In addition, selection effects could be responsible of 
increasing the scatter of the observed scaling relations 
of early-type galaxies. Early-type galaxies are usually se- 
lected by their photometric properties, being the final sam- 



ples contaminated by galaxies with di scs (e.g. SOs) which 
are not pure spheroidal systems (see IBernardi et al. 2010b . 
Color is also one of the most used selection criterion. 
However, it has been recently observed that the red se- 
quence containts a large fraction (« 50%) of red spiral 
galaxies (see, e.g., [E msellem et al. 2007 , Krajnovic et al. 2008) 
Ma hajan & Raychaudh ury 2009} IMasters et al. 20101 Sanchez- 
Almeida et al. 2011). The above contamination may explain 
differences in the reported scaling relations, since d iscs and 
spheroids follow different trends (see lShen et al. 20031) . 

In the present paper, we revisit the scaling rela- 
tions of red- sequence elliptical galaxies using a new ap- 
proach applied to one of the largest dataset of galax- 
ies available today: the spectroscopic sample of the Sloan 
Digital Sky Survey Data Release 7 (SDSS-DR7). The 
new approach of the present work is probabilistic, based 
on the spectroscopic and morphological classifications re- 
cently developed by our gro up (Sanchez Almeida et al. 20101 
Huertas-C ompan y et al. 201 1 ). This double classification pro- 
vides a robust sample of early-type galaxies avoiding the 
contamination by other morphological types. Thus, early-type 
galaxies are formed by three different classes of objects accord- 
ing to their spectra, and we study the scaling relations of these 
three populations. The differences observed are interpreted in 
terms of different formation/evolution mechanisms. 

The paper is organized as follows. We describe the sample 
selection in section^ Scaling relations such as color-mass, size- 
mass, Faber- Jackson, and the fundamental plane are shown in 
Section 3. The discussion and conclusions are given in Sections 
4 and 5, respectively. The cosmological parameters adopted 
through this paper were £l m = 0.3, = 0.7, and Hq = 70 
km s" 1 Mpc -1 . This paper is focused on elliptical galaxies (ex- 
cluding SOs). 

2. Sample selection 

Our start ing point for the samp le selection is all galaxies in 
SDSS DR7 ( [Abazajian et al. 2009[ ) with spectra of good qual- 
ity (not too close to the edges of the survey, not saturated, and 
properly deblended), apparen t magnitude in r brighter than 17.8, 
and redshift below 0.25 (see lSanchez Almeida et al. 20 101 for a 
detailed description of the sample selection). As result, the final 
number of galaxies used in the present work is 698420. 

Traditionally, early-type galaxies have been visually 
selected. This can be easily done for small and nearby 
galaxy samples. Nevertheless, visual classifications of sam- 
ples as large as the one we use are very time consuming, 
and only by involving a large number of classifiers it can 
be done in a reasonable amount of time (e.g . Galaxy Zoo 
project ILintott et al. 20081 ILintott et al. 2011b . Automated 
galaxy classification techniques have been developed dur- 
ing the last decades in order to classify large databases. 
The baseline of these techniques is that relative simple 
criteria using global galaxy properties can separate dif- 
ferent galaxy types. Thus, early-type galaxies have been 
selected accordin g to: measurements of the surface bright- 
ness profiles ([Aguerri & Trujillo 2002|Aguerri et al. 2004 



Trujillo & Aguerri 2004[~ [Gu tierrez et al. 2004 

Mendez-Abreu et al. 2008b . colors (Strateva et al. 200]}, light 
concentration and asymmetry ( Abraham et al. 1996), or some 
combination of photometric and spectroscop ic information 
(IBernardi et al. 20031 [Nigoche-Netro et al. 2010] ). 

These classifications assume that the different types of galax- 
ies are well separated in the phase space defined by the consid- 
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ered parameters. Nevertheless, this is not the case. There is a 
continuous transition between the different galaxy types. This 
makes it difficult to define a pure elliptical galaxy sample since 
contaminations from other galaxy types are always present. In 
our present work, the selection of elliptical galaxies has been 
done implicitly using a probability based approach not explored 
before. We use the full sample weighting the derived relation- 
ships with the probability of being elliptical, both morphologi- 
cally and spectroscopically. Thus even if spirals and ellipticals 
with uncertain classification are included, they do not contribute 
to the results. 

The morphological classification is based on support vec- 
tor machines (SVM). This is a machine learning algorithm 
that tries to find the optimal boundary between several clouds 
of points in an N-dimension al space (for more details, see 
Huertas-Company et al. 2008] ). The galaxies of the SDSS DR7 
spectroscopic sample were classified in four different morpho- 
logical types (E, SO, Sab, and Scd). We associated to each galaxy 
of the samp le a probability of being in the se four morphological 
classes (see Huertas-Company et al. 201 1 ). We have also classi- 
fied the galaxy spectra using the k-means clustering algorithm. 
This algorithm assumes that each galaxy spectra is a vector in a 
multi-dimensional space, with as many dimensions as the num- 
ber of spectral wavelengths. It is assumed that the vectors are 
clustered around a number of cluster centers. The k-means al- 
gorithm finds the number of clusters, the cluster centers, and 
assigns to each galaxy one of the clusters. The galaxies of the 
sample were finally classified in 17 major ASK spectral classes, 
and each galaxy has also a probab ility of belonging to all classes 
(see lSanchez Almeida et al. 2010b . 

We have combined the available spectroscopic and morpho- 
logical information by assigning a weight (Pi(E)) to each galaxy. 
This weight is the product of the probability of the i-th galaxy 
to be classified morphologically as elliptical (P/, m orph(E)) and 
the probability of the galaxy to belong to the j-th spectral class 
(^■,ask(A ; )). Note that Pi(E) would be a probability in the case 
that the two probabilities were independent. Otherwise, Pi(E) is 
just a weight which combines the spectroscopic and morpho- 
logical information for each galaxy. In the following, we will 
consider that, for each galaxy, the larger the value of Pi(E), the 
higher the probability of being elliptical. 

Early-type galaxy samples usually include both elliptical 
and SO galaxies. Nevertheless, these two kinds of galaxies 
are dynamically different. In particular, SO galaxies present 
a disc structure absent in ellipticals. This different dynamics 
can affect to the scaling relations of the galaxies. Then, dif- 
ferences observed in the scaling relations of early-type galax- 
ies could be related to different dynamics, rather than due to 
different origins and/or evolutions (see e.g. IPersic et al. 1 996; 



Bell & de Jong 2001[IShen et al. 2003UCourteau et al. 2007b . In 
order to have a sample of pure spheroidal systems, we have not 
included in Pi(E) the probability of each galaxy to be classified 
as SO. 

Galaxies classified morphologically as ellipticals do not be- 
long to all ASK spectral classes. They are basically concentrated 
in the ASK classes called ASK0, ASK2, and ASK3 (see Sanchez 
Almeida et al. 2011). The concentration of elliptical galaxies 
in only three spectral classes is just because blue galaxies with 
elliptical morp hology are very rare in the local universe (less 
than 10%, see [Kannappan et al. 2009, Sanchez Almeida et al. 



account this approximation, the weight assigned to each galaxy 
of the catalog is given by 

Pi(E) * Pi,morph(E) X (P;,ASK(A ) + P/,ASk(A 2 ) + ^z,ASk(A?)), (1) 

where P/,ask(Ao), ^,ask(Ai), and P;,ask(A 3 ) are the probabil- 
ities of belonging to ASK0, ASK2, and ASK3 spectroscopic 
classes, respectively. 

In the present work, we have assigned three different weights 
to the family of elliptical galaxies corresponding to their three 
different main spectral types (ASK0, ASK2, ASK3). Hereafter 
we will call them EASK0, EASK2, and EASK3, respectively. 
The values of the weights are given by, 



P/(EASK0) = P U morph{E) X Pi,ASK(A ), 
P;(EASK2) = Pi,morph(E) x Pf, ask(A 2 ), 
P/(EASK3) = Pi,morph(E) x P itA s K (A 3 ). 



(2) 
(3) 



(4) 



It is obvious from previous equations that: Pi(E) = 
PKEASKO) +P,(EASK2) +P i (EASK3). In all relations analyzed 
throughout the paper, all galaxies contribute to each relation with 
the above weights. Figure[T] shows snapshots of the four galaxies 
with the highest ^(EASKO), Pi(EASK2), and P,(EASK3) val- 
ues. 

In nature, we expect a gradual transition between the differ- 
ent spectroscopic and morphological galaxy types. The approach 
presented here is therefore adapted to this characteristics in the 
sense that each galaxy of the catalog belongs to the EASK0, 
EASK2, and EASK3 classes in proportion to its weight. This 
has the advantage that we have not introduced additional biases 
in the selection of the elliptical galaxy sample. Our set of galax- 
ies has the same biases as the SDSS DR7 spectroscopic catalog. 

Figure [2 shows the mean spectra of the EASK0, EASK2, 
and EASK3 galaxies. These template spectra were obtained by 
combining the ASK0, ASK2 and ASK3 spectra weighted as 
given in equations (2), (3) and (4), respectively. Notice that 
the three spectral classes correspond to galaxies showing typ- 
ical spectra of early-type systems, i.e., well defined 4000 A 
break, strong absorption Balmer lines, and no strong emission 
lines. Weak emission lines can be observed at Ha, [Nil], [SII] 
and [Oil] wavelengths. Thus, galaxies from ASK0, ASK2 and 
AS K3 classes show Ha equivale nt widths EW(Ha) < 1 A 
(see I Sanchez Almeida et al. 2010b . These weak emission lines 
in early-type galaxies have been observed before. Thus, up to 
60 % of the early-type galaxies show EW(Ha) ^ 0.5 A (see, 
e.g.,|Phillips et al. 1986|). Notice also that the emission of ASK0, 



2011). Although, they become more important in the high red- 
shift Universe (see |Huertas-Company et al. 2010| ). Taking into 



ASK2 and ASK3 galaxies is in the region of the BPT dia- 
gram ( B ald wm" et al 1 9 8 Tt traditionally attributed to AGNs and 
LINERs (Sanchez Almeida et al. 20101) . However, this region of 
the diagram is also characteristics of ioni zation provided by old 
stellar populations (IBinette et al. 19941 IStasinska et al. 2008b . 
with the EW of Ha being the paramete r to distinguish the two 
possibilities dCid F ernande s et al. 201 lb . Thus, their weak emis- 
sion lines, together with their position on the BPT diagram, 
seems to be consistent with EASK galaxies being retired galax- 
ies where the interstellar medium is ionized by hot low mass 
evolved stars. 

The ASK classification does not give a class of early-type 
galaxies without emis sion lines similar to t hose galaxies already 
observed (see e.g. IBressan et al. 2006] lAnnibali et al. 20101 
IPanuzzo et al. 201 fb . Since most of the early-type galaxies show 
weak emission lines (see e.g. |Phillips et al. 1986), the presence 
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SDSS J092553162+492013.5 + SDSS J160812I84+273416.1 + SDSS J091925. 25 + 010335. 2 + SDSS J230435.43-090906.0 



Fig. 1. Cutout images of EASKO (top row), EASK2 (middle row), and EASK3 (bottom row) galaxies with the highest weights. 



of emission lines is a feature bypassed by the automatic classifi- 
cation algorithm. In order to separate them out, one would need 
to refine the classes by re-running the classification upon spectra 
of individual ASK classes (e.g., Morales-Luis et al. 2011), but 
such work has not been carried out yet. 



Padi lla & Strauss 20081) . In our approximation, the b/a distribu- 
tion of the elliptical galaxies is given by, 



Nib I a) = 



A(b/a) 



Pm n( 



(b/a)t - b/a 
Mb/a) 



), 



(5) 



2.1. Comparison with other samples of early-type galaxies 

In this section we quantify the contamination by spirals of our 
galaxy sample. We compare the apparent axis ratios distribution 
and the luminosity functions of our elliptical galaxies with other 
samples of visually classified ellipticals. In particular, we use as 
reference the recently published eyeball classification of 14034 
galaxies by Nair & Abraham (2010). 



2.1 .1 . Apparent axis ratio distribution 

The intrinsic structure of galaxies is related to their apparent 
axial ratio (b/a). Elliptical and disc galaxies show very differ- 
ent axial ratio distribution. In particular, the axial ratio distribu- 
tion of ellipticals peaks dXb/a « 0.8, and decreases to zero at 
b/a « 0.2. In contrast, discs show an almost flat b/a distribution 



where the sum is over all the galaxies in the sample, A(b/a) rep- 
resents the bin size of the histogram, (b/a)i is the aparent axial 
ratio of the i-th galaxy, and II is the rectangle function given by, 



U(x) 



1 if |*| < 1/2, 
elsewhere. 



(6) 



in the range b/a = 0.2 - 0.8 (e.g., Tremblay & Merritt 1995 



Note that the histogram defined in Eq.(6) has an integral 
equal to 2/ ^/(E). 

Figure [3] shows this axial ratio distribution. We have also 
overplotted the axial ratio distributions of those galaxies mor- 
phologically classified as E, E+S0, and E+SO+Sab by Nair & 
Abraham (2010). Notice that the b/a distributions of E, E+S0 
and E+SO+Sab all peak at about b/a ~ 0.8. Nevertheless, the 
inclusion of discs changes the tail of the distribution. Thus, 
discs increase the fraction of galaxies with b/a < 0.5. In Nair 
& Abraham sample, 2%, 4%, and 43% of the E, E+S0, and 
E+SO+Sab galaxies show b/a < 0.5. In our selected E sample, 
only 1% of the galaxies have b/a < 0.5. This indicates that the 
tail of the b/a distribution of our E sample is closer to the tail 
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1.5 



EASK 
EASK 2 
EASK 3 




0.4 



0.5 0.6 0.7 

A [jLtm] 



0.8 



0.9 



Fig. 2. Template spectra of the EASKO (red), EASK2 (green), and EASK3 (blue) spectroscopic classes. The small panels zoom into 
of the spectral ranges of the 4000 A break (left), Mg absorption triplet (center), and Ha (right). The fluxes of the three spectra are 
in dimensionless units. They were normalized to the average flux in the g-filter bandpass. 



of the axial ratio distribution of the visually classified ellipticals. 
Then, we can conclude that our automated galaxy classification 
is not significantly contaminated by discs, which would appear 
as an extended tail towards small b/a. 

We have also overplotted in Fig J3] the axial ratio distribution 
of the early-type galaxies selected in Strateva et al. (2001) and 
Bernardi et al. (2003). As it is revealed by the tail of the distri- 
butions, both samples are contaminated by disc galaxies, being 
the color se lected sample ( Stratev a et al. 20011) the most contam- 
inated (see Ber nardi et al. 20101 for a discussion about biases in 
selection of early-type galaxies). 

The b/a distributions of EASKO, EASK2, and EASK3 
galaxies can be computed by substituting in eq. (5) Pi(E) with 
P i (EASK0),P i (EASK2), and P/(EASK3), respectively. These 
b/a distributions are shown in Fig. |4] We do not find system- 
atic differences between the different spectral classes of galaxies 
and the global elliptical galaxy sample. This fact also indicates 
that if there is some disc contamination left, it is the same for the 
three classes of elliptical galaxies defined in the present work. 

2.1.2. Luminosity function 

The galaxy luminosity function, O(L), is one of the fundamen- 
tal statistics in galaxy studies. It measures the number of galax- 
ies in a given volume per bin of luminosity. In magnitude lim- 



ited samples, it computations requires considering the fact that 
brighter galaxies can be seen further away. This is done by com- 
puting V max , the maximum volume within which a galaxy of a 
given magnitude could have been detected. Each galaxy is then 
weighted by l/V max - Thus, the number of elliptical galaxies per 
unit magnitude and volume is given by, 



O(M) 



AM ^ 



fi(E) Mi 
V 



M 



AM 



-), 



(7) 



where the sum comprises all the galaxies of the sample, AM rep- 
resents magnitude bin size, and V^ max is the maximum volume 
of the i - th galaxy. In a similar way, we can compute the lu- 
minosity functions of EASKO, EASK2, and EASK3 galaxies by 
replacing P f (E) with ^(ASKO), P/(ASK2), or P f (ASK3). 

Figure[5]shows the luminosity functions (LF) of our elliptical 
selected galaxies, and those from EASKO, EASK2, and EASK3 
galaxies. We have also overplotted the LF of those galaxies vi- 
sually selected as E, E+S0, and E+SO+Sa by Nair & Abraham 
(2010). 

Notice that the LF of our EASKO, EASK2, EASK3, all com- 
bined, follows the LF of the visually identified ellipticals (top left 
panel of Fig. 0. In addition, Fig. [5] shows that EASK2 galaxies 
dominate the LF of our ellipticals at all absolute magnitudes. The 
LF of EASK3 galaxies are close to the LF of visually selected E 
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1.0 f 



1.2 



Fig. 3. Axial ratio distribution of elliptical galaxies as com- 
puted in this work (black full line). The axial ratio distributions 
of galaxies visually classified by Nair & Abraham (2010) as 
Elliptical (red line), Elliptical + SO (blue line), and Elliptical + 
SO+Sab (green line) are also overplotted. The dashed and dot- 
ted lines represent the axial ratios of early-type galaxies from 
Bernardi et al. (2003) and Strateva et al. (2001), respectively. 
All histograms have been normalized to the same area. 



0.8 - 



0.6 - 



0.4 ■ 



0.2 I 




9 10 11 12 

log(M star ) 

Fig. 6. Color-mass diagram of the galaxies from SDSS-DR7 
spectroscopic sample (grey scale). Colored regions show the lo- 
cation of galaxies with P;(EASK0) > 0.5 (red), P,-(EASK2) > 
0.5 (green), and P;(EASK3) > 0.5 (blue). The red, green, blue 
and black solid lines represent linear fits to the color-mass re- 
lations weighted with P/(EASK0), Pi(EASK2), Pi(EASK3), and 
Pi(E), respectively. 




0.0 
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Fig. 4. Axial ratio distributions of EASK0 (red), EASK2(green), 
and EASK3 (blue) galaxies. All the histograms are normalized 
to their peaks. 



galaxies at the low magnitude regime M r < -18 (see also Sec. 
4). 

These facts also support that our rendering of elliptical 
galaxies is not contaminated by discs. 



3. Scaling relations of elliptical galaxies 

In this section we present the scaling properties of the galax- 
ies EASK0, EASK2, and EASK3. The relations were derived 
including all galaxies but weighted according to /^(EASKO), 
P;(EASK2) and ^(EASK3), respectively (see Sect. 2). 



3.1. Color-mass diagram 

Figure [6] shows the color-mass diagram of all the galaxies in 
the SDSS-DR7 spectroscopic sample. The m g - m r color of 
the galaxies was obtained from the petrosian g and r magni- 
tudes downloaded from the SDSS DR7 database. These mag- 
nitudes were amendded to account for Galactic dust attenua- 
tion and k-correction. Galactic absorptions and k-corrections 
were computed from the MPA-JHU DR7 spectrum measure- 
mentfl by subtracting the fiber magnitude and the derredening 
and de-deredshifted magnitudes of each galaxy (plug-mag and 
korjnag). Stellar masses were also obtained from the MPA- 
JHU DR7 release. They were computed following the strategy 
in Kaufmann et al. (2003). We have also overplotted in Fig. 
the regions occupied by galaxies with /^(EASKO) > 0.5, 
P/(EASK2) > 0.5, and ^(EASK3) > 0.5. These areas enclose 
the galaxies representative of the three spectral classes. Note that 
these galaxies populate the red-sequence of the color-mass dia- 
gram. 

It can also be seen in Fig. that the mass range of the three 
classes of galaxies is different. Thus, galaxies with /^(EASKO) > 
0.5 span the range 10.2 < \og(M star ) < 11.8, whereas those 
with P f (EASK3) > 0.5 have 9.8 < \og(M star ) < 11.3. In con- 
trast, galaxies with ,P i (EASK2) > 0.5 occupy the full mass range 
(9.8 < \og(M star < 12). We refer to Section |H for a detailed di- 
cussion on the mass distribution of the different classes. 

We have performed linear fits m g - m r = a + b \og(M star ) 
to the color-mass relation using P,(E), ^(EASKO), P/(EASK2) 
and P/(EASK3) as weights. The best fit coefficients are given in 
Tab.[T] Notice that for a given stellar mass, galaxies from EASK3 
class are bluer on average. 
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Fig. 5. Luminosity functions of our elliptical galaxies (top left, black open circles), EASKO (top right, red open circles), EASK2 
(bottom left, green open circles), and EASK3 (bottom right, blue open circles). In all panels we represent as reference the luminosity 
functions of E (open stars), E+SO (open squares) and E+SO+Sa (open diamonds) visually classified by Nair & Abraham (2010). 
The errors of the luminosity functions were computed by bootstraping. 



3.2. Mass-size relation 

According to the virial theorem, masses and sizes of galaxies are 
well correlated. Basically, more massive (or luminous) galax- 
ies are larger than their low mass counterparts. This correla- 
tion is the so-called mass/luminosity-size relation, and it is fol- 
lowed by galaxies in all range of luminosities and morpholog- 
ical types. However, there are important differences depending 
on the morphological type. In particular, the mass-size relation 
of early-type galaxies is less scattered than the one of late-type 
syst ems, even if ther e is a strong dependence with luminosity 
(see lShen et al. 20031) . Figure [7] shows the relation between the 
effective radius (r e ) and the stellar mass of the galaxies in our 
sample. We have taken the radius containing 50% of the total 
petrosian galaxy luminosity as an estimate for r e . Figure [7] also 
shows in colors the location of the galaxies with P l (EASKO), 
P,(EASK2), and P,(EASK3) larger than 0.5. Note that the size- 
mass relation steepens for ellipticals as compared to the full set 
of galaxies, and so, as compared to disc galaxies. We have also 
performed linear fits to the size-mass relation using /^(EASKO), 
P/(EASK2), and ^(EASK3) as weights. They are shown in Fig. 
[TJ The coefficients of the best fit \og(r e ) = a + b \og(M star ) are 
listed in TabQ] Notice that we also see deviations from a pure 

1 http://www.mpa-garching.mpg.de/SDSS/DR7 



Table 1. Coefficients of linear fits to the color-mass, size-mass, 
Faber- Jackson and FP relations for EASKO, EASK2, EASK3, 
and as well as for all of them combined. 



Relation 


Galaxy class 


a 


b 


Color-Mass 


EASKO 


0.51+0.46 


0.02+0.04 




EASK2 


0.41+0.26 


0.03+0.02 




EASK3 


0.36+0.41 


0.03+0.04 




All 


0.32+0.13 


0.04+0.01 


Mass-size 


EASKO 


-5.29+0.42 


0.54+0.04 




EASK2 


-4.92+0.26 


0.50+0.02 




EASK3 


-4.25+0.43 


0.44+0.04 




All 


-4.78+0.13 


0.49+0.01 


Faber-Jackson 


EASKO 


-0.49+0.45 


0.25+0.04 




EASK2 


-0.72+0.28 


0.27+0.03 




EASK3 


-0.95+0.45 


0.29+0.04 




All 


-0.81+0.14 


0.28+0.01 


Fundamental Plane 


EASKO 


3.65+0.26 


0.65+0.05 




EASK2 


3.65+0.17 


0.64+0.04 




EASK3 


3.37+0.31 


0.58+0.06 




All 


3.57+0.09 


0.63+0.02 



straight line at the low and high mass ends reported in previous 
works (see Bernard i et al. 20101) . 
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1.2 




9 10 11 12 

log(M star ) 

Fig. 7. Effective r-band radius (r e ) as a function of stellar mass 
for galaxies in the SDSS DR7 spectroscopic catalog (grey 
scale). The colored regions show the location of galaxies with 
Pf(EASKO) > 0.5 (red), P f (EASK2) > 0.5 (green), and 
P,-(EASK3) > 0.5 (blue). The solid lines show the best linear 
weighted fits to the relations (see Fig. [6] for the color code of 
these lines). 



\2 [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1.2 




-6.0 -5.5 -5.0 -4.5 -4.0 

log(aV) 

Fig. 8. Fundamental plane of all galaxies from the SDSS-DR7 
spectroscopic catalog (grey scale). The colored regions show the 
location of galaxies with ^(EASKO) > 0.5 (red), P,-(EASK2) > 
0.5 (green), and P;(EASK3) > 0.5 (blue). The continuous lines 
show the best linear weighted fits to the relations (colors of the 
solid lines as in Fig|6]). 



Galaxies in the EASK0 and EASK2 classes show similar 
size-mass relations. In contrast, galaxies from the EASK3 class 
follow a size-mass relation with a slightly different slope. In 
particular, for a given stellar mass, galaxies from the EASK3 
class are larger. The relation obtained for all the 3 classes to- 
gether (using P(E) as weight) is consistent with the one reported 
by Hyde & Bernardi (2009) for their early-type galaxy sample 
(logfo) = (-4.79 ± 0.02) + (0.489 ± 0.002) \og{M star )\ 

3.3. The Fundamental Plane 

Galaxies are not located randomly in the space defined by 
the effective radius (r e ), velocity dispersion (<x), and mean ef- 
fective surface brightness (I e ). They follow a tight relation- 



ship called fundam ental plane (FP; Djorgovs ki & Davis 1987 
Dres sier et al. 19871) given by, 



r e oc cr a L 



(8) 



This relation is a consequence of the dynamical equilibrium 
(virial theorem) together with the regular behavior of both mass- 
luminosity ratio and structure of early-type galaxies. According 
to the virial theorem, a = 2 and J3 = 1. There are, however, 
differences between these theoretical values and those obtained 
from observations. These discrepancies are called the tilt of the 
FP, and imply that the mass-to-light ratio is a function of galaxy 
mass or luminosity (M/L oc U). 

We have adopted the approach of plotting the FP in the r e vs 
cr 2 I~ l virial plane coordinates (see e.g. IRobertson et al. 20061) . 
In this representation the tilt of the FP can be quantified by the 
relation 



(9) 



where A = 1 means an alignment of the FP with the virial expec- 
tation. 



We use the velocity dispersion at the effective radius (cr re ) as 
an estimate of the central velocity dispersion (<x). Aperture cor- 
rections were however applied to the values given in the SDSS 
catalog. These corrections depend on the galaxy type since the 
shape of the velocity dispersion radial profile depends on mor- 
phology. As a matter of fact, early-type galaxies present impor- 
tant radial gradients in their velocity dispersion while late-type 
systems show almost flat profiles. Following our weight-based 
solution, the aperture correction applied to the velocity disper- 
sion of each galaxy (o- ije ) is given by, 



(^lor X (^<E) + Pi(S0) + Pi(Sab))+ 



i,SDSS 



xPf(Scd))), 



(10) 



where ct^sdss is the velocity dispersion of the i-th galaxy 
given in the SDSS-DR7 catalog P f (E), P f (S0), P f (Sab) and 
P/(Scd) are the probabilities of each galaxy to be clas- 
sified as elliptical, SO, Sa b and Scd, respectively (see 
Huertas-Company et al. 201 i) , and cr irCor is the aperture correc- 
tion for early-type galaxies given by Jorgensen et al. (1995): 
o~i,cor = cn,sDSS x (^) 004 , with r ap = 1.5 arcsec correspond- 
ing to the radius of the SDSS fiber. When P(E), P(S0) or P(Sab) 
are large, then the computed <x is close to the one worked out 
by Jorgensen et al. (1995). No aperture correction is taken into 
account for galaxies with a large probability of being Scd. 

Figure [8] shows the distribution of all galaxies in the \og(r e ) 
vs log(<x 2 /~ l ) plane. We have also overplotted the location of 
galaxies with P/(EASK0), P/(EASK2), and />;(EASK3) greater 
than 0.5. Figure [8] also shows the linear fits \og(r e ) = a + 
blog(cr 2 I~ l ) using these probabilities as weights. The coeffi- 
cients of these relations are given in Tab[TJ Notice that all classes 
present FPs with tilt (b ^ 1). The tilt is however larger for 
EASK3 galaxies. 



The velocity dispersion of the galaxies was corrected for the spec- 
tral resolution of SDSS 
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Fig. 9. The Faber- Jackson relationship for the galaxies in the 
SDSS-DR7 spectroscopic catalog (grey scale). The colored re- 
gions show the location of galaxies with /^(EASKO) > 0.5 (red), 
P/(EASK2) > 0.5 (green), and P,(EASK3) > 0.5 (blue). The 
solid lines correspond to the best linear weighted fits to the rela- 
tions (colors of the solid lines as in Fig. [6]). 



3.4. The Faber-Jackson relation 



The Faber-Jackson relation (FJR; IFaber & Jackson 19761) is one 
of the most studied projections of the fundamental plane. It 
shows a tight correlation between mass (luminosity) and central 
velocity dispersion for early-type galaxies, i.e. more massive (lu- 
minous) galaxies have larger central velocity dispersions. Figure 
[9] shows the FJR for all galaxies in our sample. 

Like in previous figures, we also show the locations of 
galaxies with P,(E), P f (EASK0), P;(EASK2), and P/(EASK3) 
greater than 0.5 and the best weighted linear fits log(<x) = 
a + blog(M star ). Coefficients of the best fits are given in Tab. 
[T] Again, galaxies from EASK0 and EASK2 classes show simi- 
lar FJRs. However, galaxies from EASK3 show a different rela- 
tion. Thus, for a fixed stellar mass, galaxies from EASK3 show 
smaller velocity dispersion than those from EASK0 and EASK2 
classes. We have also fitted the FJR for all the elliptical sample 
(see Tab[T]). This relation is also consistent with the one reported 
by Hyde & Bernardi (2009) for their early-type galaxy sample 
(log(cr) = (-0.86 ± 0.02) + (0.286 ± 0.002) \og(M star )). 



4. Discussion 

In the previous section we have derived the scaling relations of 
the 3 spectral classes defining our elliptical sample. Before inter- 
preting them, we compute the total weight of EASK0, EASK2, 
and EASK3 galaxies. We find that 28% of the total weight 
correspond to EASK0, 49% to EASK2, and 23% to EASK3. 
EASK2 galaxies therefore represent half of the total sample. We 
also show in Fig. [TO] the mass fraction (MF) of the 3 spectro- 
morphological classes. The MF of a given class (EASKi; i=0,2, 
or 3) and stellar mass (M star ) was obtained with the following 
expression: 



MF(EASKi;M star ) = 



2j P/(EASKi) n( 



Mstar j-M & 



(ii) 
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Fig. 10. Mass fraction of EASK0 (red points), EASK2 (green 
points) and EASK3 (blue points) galaxies in each stellar mass 
bin. 



where the sum is over all the galaxies of the sample of a 
given M star , and AM is the stellar mass bin. Clearly, the 
EASK2 class dominates in a wide range of stellar masses (10.2< 
log(M^ tor ) < 11.7). EASK3 becomes dominant in the low mass- 
end (log(M star ) < 10) while EASK0 galaxies tend to be more 
abundant at high stellar masses. In the following we will there- 
fore consider the EASK2 class as the canonical class of ellipti- 
cal galaxies, and interpret the relative differences seen in the two 
other classes. 

We have also studied the star formation histories (SFH) and 
metallicities of each spectro-morphological class. This was done 
as an independent confirmation of the results based on the spec- 
tral indexes. This was done using the data from Cid-Fernandes 
et al. (2005). They obtained the SFH of galaxies in SDSS-DR7 
spectroscopic catalog by using STARLIGHT, a code which fits 
the observed spectrum of each galaxy by a combination of single 
stellar populations (SSP) spectra of different ages and metallici- 
ties. Thus, the model spectrum is given by: 



(12) 



where is the spectrum of the j-th SSP normalized at Aq, 
M Aq is the synthetic flux at the normalization wavelength, N* 
is the number of SSPs used in the fit, Xj is the population 
vector, and the symbol <g> denotes convolution. The extinction, 
yyt, assumes a galax y-like law with a single free parameter 
(ICardelli et al. 1989b . The modelled spectrum M A is convolved 
with a Gaussian distribution, G(v*,<x*), centered at velocity v* 
and with dispersion <x*. Since all SSPs are assumed to start with 
the same mass, the components of the population vector Xj rep- 
resent the fractional mass contribution of each SSP to the model 
flux at Aq. For the SDSS-DR7 galaxies, Cid-Fernandes et al. 
used 150 single stellar population templates based on Bruzual 
& Chariot (2003) models combined according to Padova 1994 
evolutionary tracks (Girardi 1996, and references therein). They 
cover a grid of 6 metallicities (from 0.005 to 2.5 times solar) 
and 25 ages (from 1 Myr to 18 Gyr). Further details are given 
in Sect. 2.1 of Asari et al. (2007). The Code searches for the 
minimum^ 2 = £^((6^ - M^)w / i) 2 , where Ox is the observed 
spectrum and w~ l its errors. STARLIGHT uses the Metropolis 
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scheme for the x 1 minimization (see ICid Fernandes et aL"2 005 
for a full description of the code). 

The population vector Xj of each galaxy in the SDSS sample 
was computed by Cid Fernandes (2010, private communication). 
We use them to characterize mean SFHs and metallicities. 

Like in previous sections, the mean SFH and metallicities 
of each spectro-morphological class of ellipticals were obtained 
by using ^(EASKO), />;(EASK2), and P,(EASK3) as weights. 
FigurJTTI shows the results as the percentage of stellar mass cre- 
ated at each time, xj, and with each metallicity, Z. 

Details are discussed below, but note that most of the stars 
were created long ago (time > 10 Gyr) and with high metallicity 
(solar or larger). 

4.1. EASK2: canonical class of elliptical galaxies 

EASK2 galaxies represent 50% of all elliptical galaxies in 
our sample. Moreover, they are well defined as spectral 
class according to the ASK classification (see Fig. 8 from 
Sanchez Almeida et al. 2010l ) in the sense that they are well 
clustered in the n-dimensional space where each spectrum is a 
point. They out number in the linear mass regime of the scaling 
relations defined by Bernardi et al. (2010). Consequently, the 
scaling relations of EASK2 are expected to be representative of 
the "elliptical class". This is confirmed in Tab. [T] showing that 
the fitting parameters of the scaling relations of EASK2 and all 
ellipticals to be very similar, and consistent with previous pub- 
lished results. 

Concerning the SFH of these galaxies, Fig. Qj] shows that 
90% of the stars were already formed 10 Gyr ago. This is the 
espected formation history of the stel lar population of a ty pical 
monolithic elliptical galaxy (see, e.g. JCarretero et al. 2007b . We 
can conclude that our classification scheme separates as EASK2 
typical ellipticals that follow the standard scaling relations and 
star formation histories. 

In addition, the class of elliptical galaxies also contains 
by two slightly different populations: EASK0, and EASK3 
classes. This separation is a continuous transition (see Fig. 
ITOl) since these two classes are not completely clustered 
in the parameter space as the EASK2 (see Fig. 8 as in 
San chez Almeida et al. 20101) . This does not mean that they are 
outliers of the main class. Galaxies from EASK0 and EASK3 
indeed represent -50% of all ellipticals. On average, EASK3 
galaxies are bluer, larger, and less massive while EASK0 are 
slightly redder and more massive (see Figs. l6l and fTUb. 

4.2. EASK0: dry major mergers origin? 

EASK0 galaxies become a significant fraction of the elliptical 
population at the high-m ass end, where the c urvature of the 
mass-size relation appears (IBernardi et al. 20101 see also Fig. [7]). 
In a more recent work, Bernardi et al. (201 1) argue that this cur- 
vature seen in the mass-size relation can be a consequence of 
major dry-mergers becoming dominant in the high-mass end of 
the mass distribution (\og(M star ) > 1 1.3). Since EASK0 galaxies 
are more abundant at this mass, they are significantly contribut- 
ing to the curvature. We speculate that EASK0 galaxies could be 
elliptical galaxies that have experienced more major dry merger 
events than the average population of ellipticals. 

The mean SFH shows that a large fraction (about 90%) of 
their stars were already in place 10 Gyr ago as for EASK2. 
Nevertheless, they have slightly higher metallicity than EASK2 
galaxies (which may also account for the redder colors seen in 




p gal (Gal Mpc" 2 ) 

Fig. 12. Cumulative distribution function of the galaxy density of 
the EASK0 (red), EASK2 (green), and EASK3 (blue) galaxies 
in the sample of Nair & Abraham (2010). 



the color-mass diagram; Fig. 0). This could imply that EASK0 
galaxies formed stars in a more efficient way, in the sense that 
they consumed a larger fraction of the original gas leading to 
larger metallicities. Indeed, the difference of the SFH between 
EASK2 and EASK0 galaxies is positive, meaning that EASK0 
formed their stars somewhat earlier (see FigfTTl). 

Based on the above arguments, we speculate that EASK0 
galaxies could be representative of the brightest cluster galaxies 
(BCGs). These galaxies would be located in high galaxy density 
environments (mainly centers of galaxy clusters). In these spe- 
cial regions, dry mergers are normal as gas stripping efficiently 
removes the gas co ntent of galaxies in short time-scales (see e.g., 
Qui lis et al. 20 00). If the hypotheses out to be correct, our clas- 
sification scheme is able to isolate, among the population of el- 
liptical galaxies, those which have experienced more major dry 
mergers and formed stars more efficiently on average (see FigEJ). 
Moreover, even if dominant at high masses (as pointed out by 
Bernardi et al. 201 1), they are also present at all stellar masses. 

We have investigated the environment of the EASK0, 
EASK2, and EASK3 galaxies in the sample of Nair & Abraham 
(2010). Figure [T2l shows the galaxy density (p ga i) of the EASK 
classes. As explained in Nair & Abraham (2010), the environ- 
mental overdensity was taken from Blanton et al. (2005). Thus, 
the galaxy neighbors of each galaxy of the catalog were consid- 
ered as those galaxies in the magnitude range M* ± 1 and within 
5/T 1 Mpc excluding the target galaxy. Figure [121 shows that 
EASK0 and EASK2 galaxies from Nair & Abraham (2010) live 
in higher galaxy density environments than EASK3. We have run 
a Kolmogorov-Smirnov test, showing that the cumulative galaxy 
density distribution functions of EASK0 and EASK2 galaxies 
are statistically similar. In contrast, the cumulative galaxy den- 
sity distribution function of the EASK3 galaxies is statistically 
different from those for EASK0 and EASK2. 

4.3. EASK3: wet minor mergers? 

Galaxies in the color-mass diagram are thought to evolve 
from the blue clou d to the red sequence (see, e.g., 
IKauffmann et al. 20061) . The bluer colors of EASK3 galaxies 
could possibly be caused by their recent arrival to the red se- 
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Fig. 11. Top row: Mean star formation history (left pannel) and metallicities (right pannel) for EASK0 (red), EASK2 (green) and 
EASK3 (blue). Bottom row: Differences of the SFH (left pannel) and metallicities (right pannel) with respect to EASK2 galaxies. 



quence. If true, EASK3 should be poststarbust or E+A, galax- 
ies. E+A galaxies show indeed a young stellar population to- 
gether with a lack of on-going star formation. We have tested 
that possibility by determining the average value of the equiva- 
lent width (EW) of US and [Oil] lines for EASK0, EASK2, and 
EASK3. This was done by using the EWs provided by MPA- 
JHU DR7 spectrum measurements. We obtain < EW(Hs) >= 
(-0.086 ± 0.003)A, (-0.049 ± 0.002)A, and (-0.0023 ± 0.004)A 
for EASK0, EASK2, and EASK3 galaxies respectively. On the 
other hand the < EW([OII]) >= (-2.99 ± 0.03)A, (-2.25 ± 
0.02)A, and (-2.46 ± 0.03)A for EASK0, EASK2, and EASK3 
galaxies respectively. These values rule out the hypothesis that 
EASK3 galaxies are dominated by E+A galaxies. This type of 
galaxies sh ow EW(HS) > 5.0A and EW([OII]) > -2.5A (see 
IGoto 20071) . 

Cold gas accretion into a galaxy halo can reach the cen- 
tral regions of the galaxy and turn on a central starburst (see 
IDekel et al. 20091) . This is expected to happen at the intersec- 
tion of filamentary structures of the Universe. Thus, this effect 
would be significant in galaxies with masses larger than « 10 12 
M (Dekel & Birnboim 2006 ). Since EASK3 galaxies are essen- 
tially in the low mass end, we conclude that EASK3 galaxies are 
not likely rejuvened by direct gas accreation, and we are left with 
the possibility that they were formed by mergers of galaxies with 
larger gas fraction than those forming EASK2 galaxies. 

According to the structure formation scenario, mergers are 
the main driver in galaxy formation. In particular, the progeni- 
tors of elliptical galaxies could be both spheroidal or disk-like 
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Fig. 13. H/3 versus [MgFe]' of the galaxies from SDSS-DR7 
(grey scale). The asterisks show the mean values of the two 
indexes for EASK0 (red), EASK2 (green), and EASK3 (blue) 
galaxies. 



systems. The dissipative component of the progenitors plays an 
important role in determining the final properties of the rem- 
nants. Thus, during the merger process of gas-rich disk galax- 
ies, part of the gas of the progenitors is driven towards the cen- 
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ter of the remnant and feeds a central starburst that induces 
a rejuvenation of the central part of the final galaxy (see e.g. 
IDi Matteo et al. 20051) . This new generation of central stars can 
influence the structural parameters of th e galaxies and produce 
changes in their scaling relations (see IRobertson et al. 2 006, 
|Hopkins et al. 2008] ). Therefore, we expect that elliptical galax- 
ies built up after wet mergers, or those with a wet last merger, 
will show younger stellar populations in their centers than those 
formed by dry mergers. We have tested the possibility that galax- 
ies from the EASK3 class are formed by mergers with a higher 
fraction of gas than those from EASK0 and EASK2 classes. In 
this case, we expect EASK3 to have younger stellar populations 
in their central regions than average. This seems to be happening 
according to the star formation histories represented in Fig. 11. 
In addition, we carried out an independent test by computing the 
mean age and metallicity of the galaxies from EASK0, EASK2, 
and EASK3 using two Lick spectral indexes: and [MgFe]'. 
The H^ index is an age-sensitive index, whereas [MgFe]' mea- 
sures metalicity. Thomas et al. (2003) defined [MgFe]' as 



with large gas fraction can also activate central starbursts and 
produce the rejuvenation of the central regions of galaxies 
(IBournaud et al. 20071 [Kaviraj et al. 2009|). 



[MgFe]' = VMgb x (0.72Fe5270 + 0.28Fe5335), 



(13) 



where Mgb, Fe5270, and Fe5335 are Lick indexes. 

[MgFe]' is a good tracer of the total metallicity of the 
stell ar populations due to its independence of age and a/Fe 
(see IThomas et al. 20 03 ), and ensures to break down the age- 
metallicity degeneracy. The Lick indexes of the galaxies were 
also obtained from MPA-JHU web page. The indexes used 
were those measured on the restframe spectrum after substrac- 
tion of all 3<x emission linefl Figure [13] shows the H/3 ver- 
sus [MgFe]' for the galaxies from SDSS-DR7 database. We 
have also overploted the mean values of these two indexes 
for EASK0, EASK2, and EASK3 galaxies. We have deter- 
mined the stellar population properties (age and metallicity) 
of each elliptical g alaxy class by using the program R-model 
(ICardiel et al. 20031) . The result indicates that EASK3 galaxies 
are about 2 Gyr younger than galaxies from EASK0 and EASK2 
classes. This is in full agreement with the SFH seen in Fig. [TT] 
for these galaxies. It is evident from this figure that EASK3 
galaxies have a more extended SFH than the canonical ellipti- 
cals (EASK2). While only about 10% of the stars from EASK0 
and EASK2 galaxies are formed in the last 8 Gyr, EASK3 galax- 
ies formed up to 20%. 

Age and/or metallicity gradients have been observed in some 
early-type galaxies. These gradients, together with the aperture 
bias introduced by the 3"diameter of the SDSS fibers, could po- 
tentially create age differences between our classes of ellipticals. 
In order to discard this problem, we divided our galaxies in two 
groups: r e > l'/5 and r e < l'/5. In both groups EASK3 galaxies 
are about 2 Gyr younger than EASK0 and EASK2, but the sec- 
ond group is not sensitive to aperture bias since the SDSS fiber 
covers most of galaxy. Therefore, we discard the difference of 
age being due to aperture bias. 

Early-type galaxies can be produced by major mergers 
of equal-mass disks or spheroidal galaxies. Nevertheless, sev- 
eral minor mergers can also produce galaxies with pho- 
tometrical and dynamical properties typical of early-type 
galaxies. In particular, repeated minor mergers can grow 
spheroidal galaxy components with Sersic shape parameter 
n from 1 to 4 ([Aguerri et al. 200T] lEliche-Moral et al. 20061 
IBournaud et al. 20071) o r produce galaxies w ith V/cr similar 
to early-type galaxies (IBournaud et al. 20071) . Minor mergers 



see http ://w w w. mpa- garching . mpg . de/S D S S/DR7 



The structural parameters of the final remnant depend on the 
merger type, because the mass of the progenitors is deposited in 
a different way in major and minor mergers. Thus, major merg- 
ers can easily increase the mass density in the galaxy centers. 
In contrast, minor progenitors loose a large fraction of mass in 
the outer regions of the m ain progenitor (see Aguerri et al. 200 1[ 
lEliche-Moral et al. 200"6l) . Robertson et al. (2006) and Hopkins 
et al. (2008) have investigated the scaling relations of early- 
type remnants formed by major disk- like galaxies with different 
gas fractions. They have shown that progenitors with a larger 
fraction of gas produce remnants following steeper size-mass 
relation and showing larger central velocity dispersions. They 
also show that remnants of gas rich mergers are located in a FP 
with larger tilt than those from progenitors with small gas frac- 
tions. Our elliptical galaxies from EASK3 are located in a FP 
with larger tilt than EASK2 galaxies. This could be related to 
the larger gas fraction of the progenitors of the galaxies from 
EASK3. In contrast, the size-mass and FJ relations of EASK3 
galaxies show smaller slopes. This is not in agreement with 
the gas-rich major merger scenario modeled by Robertson et al. 
(2006) and Hopkins et al. (2008). We conclude that the rejuvena- 
tion of the central regions of the EASK3 galaxies does not seem 
to be caused by gas-rich major mergers. 



We can therefore conjeture that EASK3 galaxies have expe- 
rienced more disipative minor mergers than the average ellip- 
ticals. Minor mergers increase indeed the effective radius and 
decrease the cen tral velocity dispersion of the galaxies (see 
INaab et al. 2009b . Our classification scheme can therefore pro- 
vide the typical spectrum of this kind of galaxies (see Figj2j). 
The disipative minor mergers undergone by EASK3 galaxies 
could also indicate that these galaxies are located in relatively 
less dense environments (see Fig. [12]). This assumption is also 
in agreement with the percentages of the different classes of el- 
liptical galaxies reported previously. Thus, it is well known that 
more than 50% of galaxies are located in clusters or in group 
environments (see e.g. JRamella et al. 20021). This is particularly 
true for elliptical galaxies (see e.g., Dressie r 1980b . This implies 
that more than 50% of elliptical galaxies are located in clusters 
or groups environments. We have seen in the previous sections 
that EASK0 and EASK2 galaxies are located in similar envi- 
ronments and represent about 80% of our ellipticals. In contrast, 
EASK3 galaxies are located in less dense environments are rep- 
resent about 20% of the total sample. According with this num- 
bers we can say that EASK3 galaxies could be located in the 
outskirts of galaxy clusters, or in field. In these environments 
gas stripping is not active and small galaxies can retain their gas 
content. Thus, the environment could also explain the mass di- 
chotomy observed in Fig. [TO] The role of environment in the 
age and metallicity of early-type galaxies was also observed by 
Clemens et al. (2006, 2009). Similarly to the results obtained 
here, they obtained that early-type galaxies located in field were 
younger than those located in clusters. They proposed that the 
assembling processes of elliptical galaxies are independent of 
the environment, although they are delayed in the field. This 
could be in agreement with our results if disipative minor merg- 
ers events are active until later times for EASK3 galaxies than 
for the other two classes. 
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5. Conclusions 

We have isolated a sample of elliptical galaxies in the SDSS DR7 
using an automated probabilty-based spectro-morphological 
classification. We have shown that morphologically defined el- 
lipticals are basically distributed in 3 spectral classes (EASKO, 
EASK2, EASK3). For each of these three spectral classes, we 
have studied the classical scaling relations (color-mass, mass- 
size, Faber- Jackson, and fundamental-plane), as well as star for- 
mation histories, metallicities, and mass distributions. 

The typical elliptical, understood as a galaxy which follows 
normal scaling relations of ellipticals, falls in our EASK2 class. 
This class represents ~ 50% of the whole population and essen- 
tially dominates at intermediate stellar masses. The bulk of their 
stellar content was already in place 10 Gyr ago. At the high mass 
end, EASKO galaxies become more abundant. These ellipticals 
are slightly more metal rich and seem to have formed stars in a 
more efficient way than EASK2 galaxies. The tilt of the mass- 
size relation at high masses lead us to conclude that these galax- 
ies could have experienced more dry mergers than the average. 
The low mass end is dominated by the EASK3 ellipticals. They 
are bluer, larger, and have smaller velocity dispersions at a fixed 
stellar mass. Moreover, they have a more extended stellar forma- 
tion history, i.e., 20% of their stars were formed in the last 8 Gyr 
as compared to 10% for the other classes. Minor gas rich merg- 
ers could be the main driver in the evolution of these objects if 
the rejuvenation of these galaxies comes from their central parts. 
This should be confirmed (at least fot the smallest galaxies) by 
forthcoming studies about stellar population gradients. 

The environment could be at the base of the differences ob- 
served between the elliptical galaxy classes reported in this pa- 
per. EASKO galaxies would be the brightest cluster galaxies, lo- 
cated in the centers of galaxy clusters were gas can be swept 
from galaxies in short time- scales due to gas stripping. Thus, 
dry galaxy mergers would be the main drivers of the evolution 
of galaxies located in these environments. In contrast, EASK2 
galaxies could be located in less dense environments where gas 
stripping mechanism is not active and galaxies can kept their gas 
content. 

Finally, we would like to emphasize that our classification 
can isolate the spectra of elliptical galaxies with different evo- 
lutive pathways, which could not be done with classical mass 
and/or color cuts. As a consequence, the present work can be 
used as a reference for studies at higher redshift. In forthcoming 
papers we will continue dissecting the spectro-morphological 
properties of local galaxies by extending this study to other 
galaxy classes. 
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Appendix A: Dependence of the results on the 
effective radius measurements 

We have estimated the effective radius of the galaxies as the ra- 
dius containing 50% of the total petrosian galaxy luminosity. 
This radius is not corrected for seeing, which may affect the 
scaling relations presented in this work. We have estimated this 
influence by measuring the mass-size relation of our E galax- 
ies using the effective radius given by Simard et al. (201 1). This 
different estimate was obtained by fitting the surface brightness 
distribution of the SDSS-DR7 galaxies with Sesic models (see 
Simard et al. 201 1). These models provide structural parameters 
of the galaxies after seeing correction. We have studied the influ- 
ence of the used effective radius on the size-mass, fundamental 
plane, and Faber- Jackson relations, i.e., the relations involving 

Figure IA.1I shows the mass-size relation for EASK0, 
EASK2, and EASK3 galaxies as computed in Sec. 3.2, but us- 
ing the effective radius given by Simard et al. (2011). The fits 
of the size-mass relations showed in Fig. lA.ll are given in Tab. 
IA.1I Note that the relations obtained here are consistent with 
those obtained in Sec. 3.2. In particular, for a given stellar mass, 
galaxies from the EASK3 class are also larger. 

Figure IA.2I shows the fundamental plane (FP) for EASK0, 
EASK2, and EASK3 galaxies computed using the effective ra- 
dius from Simard et al. (2011). The fits to these FPs are shown 
in Tab. IA.U Notice that the behaviour is similar to that reported 
in Sec. 3.3. In particular, all classes present FPs with tilt, which 
is largest for EASK3 galaxies. In this case, the tilts of the FPs 
are closer to the virial plane (b = 1) than those reported in Sec. 
3.3. 

Figure IA.3I and Tab. IA.1I show the Faber- Jackson relations 
of EASK0, EASK2, and EASK3 galaxies derived using the ef- 
fective radius of the galaxies from Simard et al. (2011). In this 
case, the behavior and fits to the relations are the same as those 
reported in Sec. 3.4. 

We can conclude that the results inferred from the scaling 
relationships of EASK0, EASK2, and EASK3 do not depend on 
the way the effective radius is measured. 
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Table A.l. Coefficients of linear fits to the size-mass, Faber- 
Jackson and FP relations for EASKO, EASK2, EASK3, and as 
well as for all of them combined, using the effective radii cor- 
rected for seeing by Simard et al. (201 1). 



Relation 


Galaxy class 


a 


b 


Mass-size 


EASKO 


-5.97+0.43 


0.61+0.04 




EASK2 


-5.70+0.26 


0.59+0.02 




EASK3 


-4.87+0.43 


0.51+0.04 




All 


-5.54+0.13 


0.57+0.01 


Faber-Jackson 


EASKO 


-0.48+0.45 


0.25+0.04 




EASK2 


-0.72+0.28 


0.27+0.03 




EASK3 


-0.95+0.45 


0.29+0.04 




All 


-0.81+0.14 


0.28+0.01 


Fundamental Plane 


EASKO 


4.57+0.27 


0.81+0.06 




EASK2 


4.68+0.18 


0.82+0.04 




EASK3 


4.37+0.31 


0.76+0.06 




All 


4.57+0.09 


0.81+0.01 
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Fig. A.l. Effective r-band radius from Simard et al. (2011) as a 
function of stellar mass for galaxies in the SDSS DR7 spectro- 
scopic catalog (grey scale). The colors and lines are as in Fig. 
7. 
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Fig. A.2. Fundamental plane of the galaxies from SDSS-DR7 
spectrosocpic catalog (grey color) using the effective radius 
given by Simard et al. (2011). The colors and lines are as in 
Fig. 8. 
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Fig. A.3. The Faber-Jackson relation of the galaxies from SDSS- 
DR7 spectrosocpic catalog (grey color) using the effetive radius 
given by Simard et al. (201 1). The colors and lines are as in Fig. 
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